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ABSTRACT 


A simple formula for the transition probability for electron exchange 
between unlike ions and atoms is established within the adiabatic 
approximation by employing the Linear Combination of Atomic Orbitals (LCAO) 
method. The formula also involves an adiabatic parameter, introduced by 
Massey, and thus the difficulties arising from the internal energy defect and 
the adiabatic approximation are avoided. Specific reactions 
Li^ + + H >Li ++ + H + and Be^ + + H ->-Be^ + + H + are considered as examples. The 
calculated capture cross section results of the present work are compared with 
the experimental data and with the calculations of other authors over the 

7 Q 

velocity range of 10 cm/sec to 10 cm/sec. 
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I. INTRODUCTION 

Charge transfer collisions provide an interesting three-body system which 
is amenable to theoretical analysis in addition to their importance in the 
studies of fusion plasmas, astroph.ysi cal processes, and more recent studies of 
short wavelength lasers. Traditional interest still looms as important 
applications in heavy ion accelerators and radiation-induced chemistry as well 
as charged particle transports in matter. 

In the past, the perturbed stationary state (PSS) method has been applied 
to study a symmetric resonance process such as A + + A +A + A + . Various 
theoretical approaches rest on obtaining stationary molecular states 
(refs. 1-3). These results are usually in good agreement with their 
experimental counterparts. Encouraged with this success, in the above 
resonant process, extending these calculations to the collisions between 
unlike ions and atoms is presently studied. In this context, the reaction of 
the type 


A Z + B -► A Z " q + B q + AE (1.1) 

in which AE is the internal energy defect would be of great interest. 

Unfortunately, the formulation of the problem of the non-symmetri cal 
charge transfer collisions of equation (1.1) is rather difficult using the 
procedures similar to the ones applied to the symmetric case. The correct 
procedure should involve the expansion of the wave function in terms of 
molecular wave functions. Rapp and Francis (ref. 4) have pointed out that the 
connection between the molecular and atomic approach in non -symmetrical charge 
transfer is not easily obtained. These authors in their work abandoned the 
molecular expansion, and instead, within the lowest order approximation, they 
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replaced the molecular wave functions and energies by atomic orbitals and 
atomic energies. 

In this paper, we adopt a correct molecular expansion procedure by 
employing the Linear Combination of Atomic Orbitals (LCAO) method. The non- 
symmetrical charge transfer process of equation (1.1) is calculated for low 
and intermediate energies using an adiabatic approximation together with 
Massey's (refs. 5 and 6) adiabatic criterion resulting in a limited collision 
range. A simple formula for the transition probability is obtained under the 
two-state approximation. The formula is expressed in terms of overlap 
integrals, exchange integrals, and Coulomb integrals. The internal energy 
defect AE and the collision range R c also appear in these formulas. 

Specific reactions Li 3+ + H +Li ++ + H + and Be 4+ + H +Be 3+ + H + are 
considered as examples. The results are compared with the experimental data 
and with the calculations of other authors, and a brief discussion is 
presented. All these overlap integrals, exchange intergrals, and the Coulomb 
integrals can be calculated analytically by employing the hydrogenic wave 
functions. The simplified picture presented in this paper gives insight into 
the essential nature of the charge transfer process. 

II. QUASI-CLASSICAL APPROXIMATION 

As is known, the passage from quantum mechanics to the limit of classical 
mechanics is similar to the passage from wave optics to geometric optics 
(ref. 7). The field components U in the electromagnetic wave can be written 
as U = ae 1 ^, where a is the amplitude, 4> i s the phase and i is the square root 
of minus one. The limiting case of geometric optics corresponds to a short 
wavelength i.e., to a large value of the phase <(>. Similarly, in the limiting 
case of classical mechanics, one writes the electronic wave function as 
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¥ = a e^* 1 where a is a slowly varying function, S is the action which 
takes on large values, and h is the Planck constant divided by 2ir. 

Henceforth, we will use atomic units in which h = 1. This use of the action 
variable is determined after considering the analogy of the principle of least 
action in classical mechanics and Fermat's principle in optics. As is known, 
the action function is given by 


S = - Et + E J p i dq i 


ER 


~ ~ — + Z J p i d( li 


( 2 . 1 ) 


and becomes large for small velocity u. Thus, Landau (ref. 7) stated that 
= ae^ is the wave function of an "almost classical" or "quasi -classical " 
physical system. Furthermore, by absorbing the second term of the above 
equations into the quantity a' and retaining the adiabatic approximation 
(i. e. , u extremely small), the wave function becomes 




( 2 . 2 ) 


where a' is a slowly varying function of time t. Let us now apply the above 
to the capture process of present interest. 

Consider the situation in which the projectile captures an electron from 
the target. Then is the wave function of the electron long before the 
collision and is essentially the electronic wave function of the target. On 
the other hand, after the collision the wave function of the electron would be 
the electronic wave function ip ^ °f the projectile. These states are not 
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stationary states in the sense that S - Et is an extremum along the 

trajectory. The stationary states of this system are the two molecular states 

denoted by <|> and . (In a symmetric case, g refers to a gerade state and u to 
y u 

ungerade state.) Thus, according to equation (2.2), within the adiabatic 
approximation, the wave function at anytime t can be expanded as 

♦ = A(t) * g exP(-i i! t o E g dt) + B(t) + u ex P(‘ i J-to E u dt) {2 * 3) 

where E and E u are eigen energies; i|» and <|> are the corresponding two eigen 
functions. (Two-state approximation is employed.) Equation (2.3) is similar 
to the one usually derived under the PSS method (ref. 8). 

The quantities A(t), B(t), i|>g» and are slowly varying functions of 
time. The conditions under which A(t) and B(t) can be treated as time 
independent are discussed in the appendix. Under the adiabatic approximation, 
one can ignore the terms 4< g and relative to the other terms in the 
pertinent equations. This leads to the conclusion that the quantities A(t) 
and B(t) are independent of time. 

For simplicity, consider a fully stripped ion of atomic number Z, moving 
slowly through the target hydrogen atoms. Let r^ be the distance of the 
electron to the target nucleus and the distance from the projectile nucleus 
(fig. 1). The i nternucl ear distance of the projectile and the target is 
denoted by R. The Hamiltonian of the system with the electron and two nuclei 
in atomic units is 






(2.4) 
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Wave functions of interest are 

*r77 ex P (- r i) 


(2.5) 


'l> 2 = A n Z n ' 1 r^" 1 exp(-Zr 2 /n) (2.6) 

7 3/2 2 2n " 1 

A n = Tir (2n) ! n 2n+1 (2 * 7) 


The is the wave function of the electron when it is localized near the 
first nucleus (target), ignoring the presence of the fully stripped projectile 
ion. On the other hand, is the wave function of the electron when it is 
localized near the projectile nucleus, ignoring the presence of the target 
nucl eus. 

As noted earlier, in arriving at equation (2.3), the wave function at any 
tine is given by the linear combination of two quasi -mol ecul ar stationary 
states \j>g and t u • These states are determined by the variational method by 
employing the Hamiltonian H with their form given as 


\f> = 

9 


<|> = 
Y u 


('J' i + ’J' ? ) 

7TT 2S 1 c 


- - (+1- to) 

72 + 2S 1 £ 


( 2 . 8 ) 

(2.9) 


where S = J dt is the well known overlap integral. The eigen value of 

the Hamiltonian is the potential function for the nuclear motion in the Born- 
Oppenheimer approximation. We note that the charge transfer problem is 
similar to the usual potential crossing problem. 
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The quantities A(t) and B(t) in equation (2.3) are determined from the 
initial condition <|> -»■ ^as t * -t to be satisfied by the value 


A = j ✓ 2 + 2S 
B = |/2 - 2S 


( 2 . 10 ) 


( 2 . 11 ) 


Letting t -»• «, using equation (2.3), we obtain 


* = C'{ *j+ 1>2 


( 2 . 12 ) 


where 


C 


„ - I 
2 " 2 




E n dt l 
J U J 


(2.13) 


Hence, the transition probability is given by 


where 



. 2 


sinn 


" - ill ( £ g- E u ) dt 


(2.14) 

(2.15) 


For any symmetric resonance process, the quantity aE = E - E^ plays the 
key role in determining the transition probability. We will determine it by 
the variational method. We note that the molecular wave function <f> can be 
obtained by the LCAO method as 


<t> = C^-t- C 2 \p 2 
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The expectation value of the Hamiltonian H is 

_ _ J H4> dr 
6 " ■' y 

J <Tdr 

where 

H ij = J ** H 4-jdt (2.17) 

The stationary states are obtained by the stability condition 


Hi 1 + Cl C 0 ( H, n + 


„1 2 _\ 1 2 21 
C 2 + C 2 + 2 CjCgS 


) + C« H 


2 22 


(2.16) 


3e__ _ if 

rcy “ ac 



(2.18) 


which leads to 


C 1 (H 11' e) + C 2 {H 12“ eS) = 0 (2 * 19) 
C^Hgr eS ) + C 2 (H 22 - e) = 0 (2.20) 

For the symmetric case (Z=n=l) we have H 11 = H 22 = H. Since H is a Hermitian 
operator, and are rea l functions, thus H 12 = H 21* The secular equation 
of equations (2.19) and (2.20) is 

(H - e) (H - e) = (H 12 - eS) 2 (2.21) 

which leads to two eigen energies 

_ H ' H 12 
G l' 1 - S 


( 2 . 22 ) 
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H + H 


£ o“ 


12 


2 1 + S 


(2.23) 


from which 


e 2” e 1 = 2 


(H 12 - HS) 

1 - s 2 


(2.24) 


Using the Hamiltonian and the wave functions (Z=n=l ) results In 


H = E l,0 + J 


(2.25) 


H 12= E 2,0 S + K 


(2.26) 


where E^q and E 2j q are the atomic ground state energies of the electron in 
the target and projectile orbitals, respectively. The S is the overlap 
integral given (ref. 9) as 

S = J ♦! ipj dr 

(2.27) 

= (1 + R + ^ 3 ) e' R 
and K is the exchange integral given (ref. 9) as 


S . *1*2 * 
1 


K =F-J r 


= ( ^ - | R) e" R 


(2.28) 
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and J is the coulomb integral given (ref. 9) as 


J 


. 2 

1 , *1 „ 

■ r : I 7 - d ' 


1, -2R 


= (! + p) e 


(2.29) 


Thus , 


h 12 - HS = K - JS 

so that, 

(f - 3 ) «p(-R ) - (| + 4 * S| + -^_exp(-3R) 
1 - (1 + R + -^) 2 exp (-2R) 

For, R >> 1, we have from the above: 


(2.30) 


(2.31) 


e 


2 



(2.32) 


which is on the order of the result given by Firsov, Landau, and Herring (see 
discussion in ref. 10) of -4Rexp (-R-1). We now consider the extension of 
these results for symmetric resonance transfer to non-symmetri c resonance 
transfer. 


III. NON-SYMMETRICAL RESONANCE TRANSFER 
In recent years, because of its usefulness in many areas of physics, the 
reaction A^+ B -*• A^” q + B q + AE and, in particular, the fully stripped ion 
charge transfer process have received great attention. Considerable progress 
has been made in both the experimental and theoretical studies of this 
problem. The unitarized Distorted Wave Approximation (UDWA) of ref. 11 has 
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been a successful numerical method, since it covers a wide energy region. For 
extremely low energies, the classical Barrier model of ref. 12 and the 
absorbing sphere model of ref. 13 have been of interest to experimental 
physicists because of their considerable simplicity. The other numerical 
models, on the other hand, involve complicated calculations. Still, the early 
Landau-Zener formula is widely used by many authors, since it has been studied 
fully for a long time (refs. 14-17). However, the Landau-Zener formula has 
its limitations. In this model the transition is assumed to be limited to a 
narrow region around the crossing point R C p. The potential energy separation 
varies linearly with the internuclear distance in the narrow region, namely. 


(H.- H f ) = Y (R - R„J 


CP' 


The coupling matrix element has a constant value £ . 

In 1962, Rapp and Francis (ref. 4) studied the charge exchange between 
unlike atoms. They pointed out that the connection between the molecular 
states and the atomic orbitals is not easily obtained. Therefore, as a lowest 
order approximation, they replaced the wave functions of the stationary 
molecular states by the atomic orbitals. This connection was found by 
employing the LCAO method. The physical difficulties arise when one 
simultaneously retains the internal energy defect and the infinite interaction 
range within an adiabatic approximation or a quasi -classical approximation, 
because such a situation renders transition probabilities small. Actually, 
the infinite interaction tine does not exist at all, since there are so many 
target atoms in the material that the average interatomic distance probably 
will be a suitable estimate of the upper limit of the effective interaction 
range. After such an interaction time, the projectile undergoes another 
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collision process. Massey (refs. 5 and 6) introduced an adiabatic criterion 
based on the following argument. Classically, the transition probability is 
proportional to the ratio of transition time and the collision time. Since 
now there exists an energy defect, AE, the transition time can be estimated 
as 1/A E. Meanwhile, the collision time is proportional to an adiabatic 
parameter divided by the velocity of the projectile. Only when the transition 
time can be comparable to the collision time does there exist a significant 
transition probability. The adiabatic parameter thus cuts off the interaction 
time. This parameter is of the order of the atomic dimension. Hasted pointed 
out in ref. 6, that "the analysis of a large volume of experimental data leads 
to surprisingly small probable errors (10-20 percent) with a value R c = 7A 
common to different types of reactions." This adiabatic parameter of 7A is 
the range of internuclear separation over which the charge transfer is deemed 
possible (ref. 11). Actually, the cross sections were measured as functions 
of projectile velocities. Surprisingly, the velocities corresponding to 
maximum cross sections v m were almost a linear function of aE . The slope 
leads to an adiabatic parameter R c = 7A . We apply this adiabatic parameter 
as a "cut-off" range in our calculations. 

The capture cross section is thus determined in terms of the adiabatic 
parameter R c , the energy defect AE, and the velocity of the projectile. It 
also depends on the atonic wave functions of both the projectile and the 
target and the interaction between them. In the non -symmetrical case, the two 
stationary molecular states are still denoted by and \|> u , where the 
subscripts g and u are used here only for convenience. These states assume a 
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more general form than the corresponding symmetric resonance case. These more 
general states are defined as 


♦ g* C 1 1 + f 1^2^ 

and 

i|> u = Cj {4 , j + ^2^2^ (^*2) 

where Cj and Cp are determined by the normalization condition 

1 = C 2 (J ^ dr + ?f X J f 1 f 2 dr + f j J dr} (3.3) 


so that 


C 


1 


1 

(1 + 2Sf t + f j2) 1/2 


(3.4) 


Simi larly. 



1 

(1 + 2Sf 2 + f 2 2) 1/2 


(3.5) 


The quantities S, fp and f 2 are functions of internuclear distance R. 

Under the adiabatic or quasi -classical approximation, the quantites A(t) 
and B(t) in equation (2.3) can be treated as time independent. Thus, 

equation (2.3) with equations (3.1) and (3.2) becomes 

♦ = A f 1 i|» 2 ) exp (-1 J* t E dt) 

o y 


+ B C| (* x + f 2 i> 2 ) exp (-1 j[ t E u dt) 


(3.6) 
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Since ij> must tend to when t -*•- t Q , we obtain from (3.6): 


A' 



(3.7) 


and 



(3.8) 


with A 1 = A Cj and B' = C * j 

The quantities A', B 1 , fj, and f 2 are all functions of internuclear 
distance R. Also, these quantities should be symmetric. Thus, A ' (-t Q ) = 
A'(t 0 ), B'(-t 0 ) = B'(t 0 ), fj(-t 0 ) = f j (t Q ) , and f 2 (-t Q ) = f 2 (t 0 ), since the 
variational method is invoked, by which they depend only on the distances 
between the projectile and target nuclei. Also, since the quantities A and B 
can be deemed to be time independent 


A'f 


1 




Equation (3.6) can now be written in a general form as 


(3.9) 


♦ = CJ * 1+ 


(3.10) 
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where, 



(3.11) 


where the coefficients are to he determined by the stationary phase method 
applied to the expectation value of the Hamiltonian. Hence, the transition 
probability is given by 


W 12 = 4 ( ~7~T -) 2 sir,2ri (3.12) 

L£ 2 T 1 

where 

” = |j-f 0 < E g- E u> dt < 3 - 13 > 

Notice that the transition probability of (3.12) differs from the 

f f 

1 o i O 

symmetric case by an additional factor 4 (-2 ^ — -) . This factor is 

T 2" T 1 

determined by the LCAO method, i.e., by the solutions of the basic equations 
(2.19). 


CJ (H r e) + C£(3-eS) = 0 (3.14) 

and 

CJ (3 -e S ) + C" £H - 2 e) = 0 . (3.15) 

where e = H^= J 2 i because H is Hermitian and ^ and ^ are -real; Hi = H n and 


h 2 " H 22* 
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In order to obtain a simple and explicit solution for the transition 
probability, the values of B, S, and e will be discussed. When the 
internuclear distance is not too small, the quantities -6 and S are of the 
same order, whereas e is about the order of -0.5. Thus, by replacing the 
quantity B - eS in equations (3.14) and (3.15) by B , the problem is amply 
simplified. Later in this section, we will explain in detail the 
justification for this approximation. Thus, equations (3.14) and (3.15) can 
be written as 

Cj (H r e) + B = 0 (3.16) 


and 




Co (Ho- e ) = 0 


(3.17) 


which give 


(H r e ) (H 2 - e) = (-*-B) 


(3.18) 


For the case Ho - > 0, the above equation when solved for e yields 


e = 


h 2 + h i H i ) Bi 

2 * 2 ^ * ' A J 


(i + 4) 


(3.19) 


where 


H 2- H 1 


(3.20) 


Bj= { 1 - (B/2A) - 1} A 


(3.21) 
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When the internuclear distance is not too small, 3 is rather small. Also, in 
the non -symmetrical process, except for a resonant channel, and H 2 are 
usually not very close to each other. Thus, 


H l- H 2 


(3.22) 


Hence, under the above condition, can be simplified as 


B ~ I 6 — 
B 1 4 


The equation (3.19) leads to 


e r H r B i 


z ? = H 2 + B 1 


The equation (3.16) for e = becomes 


CJ (H x - H 2 - Bj) +|C"B = 0 


thus , 


C" 

• - _ 2 _ 

r 


(H r H 2 - Bj) H r H 2 

fj2 " TIT 


For e = e 1 , the equation (3.17) becomes 


C" 

L 1 0 


+ C^(H 2 - H 1+ B x ) = 0 


(3.23) 


(3.24) 

(3.25) 


(3.26) 


(3.27) 


(3.28) 
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Thus, 


* C 2 - 6/2 - 6/2 
' Z~ C'i ~ H 2 - Hj+ Bj " H 2 - Hj 

Thus, the transition probability is 

w l 2 = 4 sin 2 n 


(3.29) 


(3.30) 


where 


A = _L 1- 

H f f 

' 1 2 


6/2 H 2" H l + B 1 

H 2" H l" T 1 V/2 


(3.31) 


and 


n » \ AT dt = \ J^> { (H 2 - Hj) + 2Bj} dt 
o o 


(3.32) 


In 


the case Hj* H^, or where Hj and H 2 are not close to the same value, 


applying (3.22), and remembering that Bj is of the order of e , we have 


H,- H.+ B. H,- H. 

A ~ 2 1 1 -. 1 1 


6/2 


6/2 


(3.33) 


and the transition probability is 

H 2 

12 2 

W 19 = — — , si r\ n (3.34) 

(H 2 - Hj)^ 

This is the required simple transition probability formula which is 
determined by H^ 2 , (H 2 - H^), and the value of n. The quantity Hj 2 is closely 
related to the overlap integral. Physically, it is reasonable to expect 
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higher transition probability with more overlap. The quantity H^- Hj is 
strongly related to the energy defect AE (at a large distance). Obviously 
there will be less transition probability with more energy defect. 

For the case Hj = H 2 , we have from equation (3.18) 

(H - e ) 2 = (S/2) 2 (3.35) 

which yields 

e = H ± 6/2 (3.36) 

f 1= 1, f 2 = -1 and A = 2 (3.37) 

2 

The transition probability for this case is given by W = sin o» which is the 
familiar formula (2.14) of the symmetric case obtained in Section II. 

Let us now discuss the assumption of the replacement of the quantity 
6 - Se by 6/2 in equations (3.14) and (3.15). Initially, this approximation 
seems somewhat arbitrary. However, 6/2 is not strictly equal to 6 - Se in 
equations (3.14) and (3.15). Recall that the equation (2.19) is derived from 
the variational principle, which implies that the LCAO molecular state 
function is made the most stable state at that particular interatomic 
separation. On the other hand, equations (3.14) and (3.15) were altered by 
the replacement 6 - Se 6/2, which implies that a set of molecular states 
are chosen that, although not the most stable, nevertheless approach the most 
stable states for large separation. In other words, a set of quasi -stati onary 
molecular states are chosen. In this way we may refer to a quasi -stationary 
molecular state expansion procedure. Actually, in chemistry, at a large 
distance, the molecular states are not stable at all. However, we are more 
interested in looking for a set of mathematically permitted basis states. 
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Probably, this is one of the reasons why the atomic expansions are widely used 
(refs. 1, 2, 4, and 18). Thus, it is more important that these molecular 
states be orthogonal. Fortunately, our quasi-stationary molecular states 
approximately satisfy the requirement of orthogonality. In the following, we 
will explicitly show how the orthogonality requirement is approximately 


satisfied. Consider 


*q = W f l+2) 


* u = C 1 (*1 + ^2^2^ 


(3.38) 


where 


(1 + Sf 1+ fp 


i/2 and c i = 


(1 + Sf 2 + i 2 ) 


c h u u r !• i 

C 2 ^ H r H 2 C 2 

V z*- - ~or and f 2 Trp 


H 2“ H 1 


/ ip ip dr 
J g u 


= CjCj J (♦jH- f 1 'f» 2 ) (♦•!+ f 2^ 2^ dT 


Cl Ci (J + J (p + ^2^ *1*2 dT ^ 


Remembering fp 2 = - 1, we have 


(3.39) 


J * g = c i c i J ( f ! + f 2 ) * 1*2 dT 


(3.40) 


For the case *H 2 , when the internuclear distance R is not too small 

usually |H 2 - Hjj >> 1 0. 5 p|; thus jfj >> 1 and | f 2 1 << 1. Also, fj + 


f 2 * f l 


" TT 72 + f and 

(1 + 2Sf J + fp 1 


2 1/2 

1 (1 + 2Sf + f ) 

2 2 
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One sees from these that 


T i 

J 1^* 'fj- J * 1 t 2 <lT = i '*'l^2 < ^ T = S 


Now, the overlap integral of the two atomic orbitals between the 
projectile and the target is small except when the internuclear distance 
approaches zero. Thus, J ^ u dr - 0. 

For the case Hj = H 2 , we have fj = 1 and f 2 = -1. Thus, 

1 *8*3* ■ — - ^T/7 1 ( f l* f 2 ) *1*2* * 0 


For the case H,,, we have f j « 1 and f^ “ -1. Thus, 


j ♦ g^A - (fj + fg) S - 0 

We observe that for the cases Hj = H 2 and H^, because of the 
cancellation between fj and f 2 , the orthogonality is preserved. From the 
orthogonality point of view, this expansion is at least not worse than the 
atomic expansion. 

The capture cross section o, according to the impact parameter method, is 
given by 


a = 2ir J P| bdb (3.41) 

*"0 

(see fig. 2) where b is the impact parameter, and P refers to the transition 
probability. Usually, P takes the value P^ »; but as discussed earlier, an 
adiabatic parameter will be employed to evaluate the cross section. The 
transition probability P is given by equation (3.34) as 
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H 2 

P = — — ? sin 2 n (3.42) 

(h 2 - h 1 )^ 

where all the quantities appearing in this equation are functions of the 

? I? 1/2 

internuclear distance R = (b + ^ R ) ' in which R c is the adiabatic 

parameter defined by 2 o t = R c « The value of the parameter R c is chosen to 
be 7A. 


IV. CHARGE TRANSFER IN Li 3+ + H 

Let us evaluate the magnitude of the transition probability for the 
following reacton: 

Li 3+ + H -»> Li 2+ (n = 3) + H + 

The magnitude of the transition probability is 


P = 
o 


H 12 

W 2 - HjT 


4 (|) 2 


(Hf H 2 ) 


(H r H 2 ) 


(4.1) 


(Hj= H 2 ) 


where 


A -0.56 + ill2- H l + _ 8 l) 

A " TH'- H 2 - Bj + O'. 56 "" 


H 12 = j I# 2 dr 

•I 

- E 2,„ s + f- K 


77 + f > *2 * 


(4.2) 


I 
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H r 1 *i<- -F^ - rf + #> fi * 


E l,0" 1 J l + R 


" 2-1 


" E 2,n ' J 2 + R 


(4.3) 


(4.4) 


The quantities 3, H^, and n are calculated by evaluating the overlap 
integral S, the exchange integral K, and the Coulomb integral J, defined as 
follows: 

M? dr 

S = J M 2 dr, K = J , and 


1 


2 2 

*f fi 

J i = J 7 - dr , J ? = J -r- dt 

1 r 2 < r i 


The integrals S and K can be evaluated by introducing well known 
transformations of ref. 9 (fig. 3): 


^ = R ^ 2"*" 1 ) • ^ ^ ^ ^ 


(4.5) 


and 


v* ~ p ( r 2 “ ^ j ) > “1^ y ^ 1 


dr = ~ (X 2 - y 2 ) d ♦. 


(4.6) 

(4.7) 
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For our special case, Z 2 = 3, = 1, and n = 3, we have 


S = J ♦j+2 


dr 


= 18 ( 1 t) 1/2 Ai|) 3 (|) 2 i (X + u) 2 (X 2 -u 2 ) e _RA dXdu 


3 ^ 2 ; '■ 2 


Now, from table 1, the value of Ag - 0*013; thus, we obtain 


S = 0.013 (1.6 R 4 + 8R 3 + 24 R 2 + 48R + 48) e" R 


(4.8) 


Similarly, the exchange integral is 


K = 0.026 (4R 3 + 8R 2 + 12R + 12) e" R 


(4.9) 


and 


8 - I ' , 2 <h 


SZ f l*2 

" E 2,nS + % " J --- dT 


= -0.RS + - K 


The calculations of the Coulomb integrals and are done in the 
following way. (See fig. 4.) 


J 


1 


i c ? 

_1 dj. _ 2n_ , -2r sin 9 d9 r dr 

r 2 11 ( R 2 + r 2 - 2rR cos e) 1 ^ 2 
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2 J e _2fc)R p 2 R 2 cfc> re_^_ L . j_ + P. ~ >] 
J w P p J 


(4.10) 


where p = r/R. Finally, one obtains 


i ad 2 / i 1 2. i® -2fc>R . , 

J l= 4R { J n e pcfc+Jje^pdp) 


1 /, A 1, -2R 

- R - (i + R ) e 


(4.11) 


One notes that 0j»- (R > l)and 1 (R -► 0). The quantity J 2 is calculated 
as follows (fig. 5): 


J 2 = i T- dT = * ( A 3 x 9 > J e 


-2r si n 0 d9 r dr 

( R 2 + r 2 - 2rR cos ©) 


1/2 


(4.12) 


= ^ - (0.043 R 4 + 0.215R 3 + 0.645R 2 + 1.29R + 1.61 + e‘ 2R 

For the calculation of n, one notes for the Z = 3 and n = 3 case 

(resonance) that „ 

Hj= -0.5 - 3J X + 


H 2 = -0.5 - J 2 + | 


B l" ^1 + (-J-) 2 -lJ * A ’ 24 ( H 2“ h i) 
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Thus, 


n = \ J*£ AE (R 1 ) dt = j } t o t (H 2 - Hj+ 2B } ) 


dt 


(4.13) 


where, R 1 = „2 + ^2 


The capture cross section, o, is given by 


o = j“ P(R ) 2ir b db (4.14) 

where P = P Q (R) sin^n in which P Q is taken from equation (4.1) and 

p p 1 p 

H j 2 » h 2“ h 1* and B 1 are a ^ functions of R c = b + j R c . 

V. RESULTS AND DISCUSSIONS 

The results of the calculations for the reactions Li +++ + H -*-Li ++ + H + 
are exhibited in figure 6. They are compared with the results based on the 
Unitarized Distorted-Wave Approximations of Ryufuku and Watanabe (ref. 18), 
the ten molecular-state calculations of Kimura and Torson (ref. 19), and the 
calculations of Stollberg and Hai-Woong Lee (ref. 17) based on the Landau- 
Zener formula. The results are also compared with the experimental data of 
Seim et al. (ref. 20). Also, the results of the calculations for the 
reaction Be 4+ + H + 8e^ + + H + are exhibited in figure 7 and are compared with 
the UDWA approximation and Exponential Distorted-Wave Approximation 
calculations of Suzuki et al. (ref. 21) and with the molecular orbital close 
coupling calculation of Harel and Salin (ref. 22). No experimental data were 
available for this process. 
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In figure 6, both the results of our calculations as well as those based 

on the Landau-Zener formula do not appear to show the tendency of the cross 

section to increase toward the high-velocity region. This is partly due to 

the implied employment by us of the Born-Oppenheimer approximation in which 

the effect of the nuclear kinetic energy terms are ignored in the 

8 

Hamiltonian. For a large projectile velocity (u^> 10 cm/sec), the Born 

Oppenheimer approximation should break down. When nuclear motion is 

i ncorporated, there will appear an additional electron kinetic energy term 

because of the relative motion between the projectile and the target. This 

becomes more obvious when our coordinate system is set on the projectile. 

1 ? 

Then the kinetic energy of the electron becomes ni e ( u 2 + * w ^ ere u 0 the 

velocity of the electron in the laboratory frame. Hence, the average kinetic 
1 2 2 

energy becomes y m e ( u 0 + u 2^* A significant additional kinetic energy term is 
obtained when is comparable to the electron velocity that is the first Bohr 

velocity. Also, if we choose the laboratory frame, the electron will receive 
an additional velocity from the momentum transfer. Thus, because of the above 
reasons, some higher states will be excited, causing a significant 
contribution to the total cross section. 

Let us now discuss the general energy dependence of the cross section. 

One expects that toward the high velocity region, the cross section should 
decrease as a result of the increasing relative velocity between the 
projectile and the target. In the extremely low energy region, the cross 
section for these processes should not rapidly approach zero with a decrease 
in u, but rather approach a not too small constant. The calculated cross 
section in the energy region of interest (below a few keV) approaches a 
constant value due to the term H 2 - Hj in the argument of the sine function. 
This makes the value of the sine function very large, and the transition 
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probability oscillates rapidly. The summation over the impact parameter 

2 

approaches an average number, sin n - 0.5. When u becomes larger, there will 
be a slow oscillating behavior, but this velocity region is beyond our 
adiabatic approximation. 

In the symmetric case, since Hj= only when the impact parameter b is 
small does the transition probability exhibit the rapidly oscillating behavior 
different from the non-symmetric situation. 

Finally, experimental physicists (refs. 23 and 24) have recently noticed 
that below a few keV, the capture cross section shows velocity independent 
behavior in helium and lithium materials. It is hoped that our results will 
help in classifying these and other relevant capture cross section 
observations. 


CONCLUDING REMARKS 

A semiclassical approach is used to derive a quasi -stationary state 
method to approximate the charge transfer cross section. The result is only 
slightly more complicated than the usual Landau-Zeuer formalism and greatly 
simpler than the coupled channel and various distorted wave formalisms. 
Unlike the Landau-Zener result, we predict nearly constant cross sections 
below a few KeV in agreement with recent experiments and the more complex 


formalisms 
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APPENDIX 


In this appendix, the conditions under which the quantities A(t) and B(t) 

can be treated as time independent in the equation (2.3) are obtained. The 

central point in this treatment is the observation that when u -*■ 0, compared 

iS 

with the phase change, a is a slowly varying function In f ~ a e . As 
discussed in Section II, this leads to 

* = A (t ) * g e _i ^ Eg dt + B(t ) t e’ 1 ^ Eu dt (B.l) 

We have shown previously that ¥ and are slowly varying functions and 
had assumed that A and B vary so slowly that they can be treated as time 
independent. Now, let us estimate how slowly and ¥ u are varying and under 

what conditions A and B can be treated as time independent. Having done this, 
we will know the conditions under which our analytical model of this paper 
works. 

The wave functions of the stationary states can be formed by linear 
combination of atomic wave functions, i.e.. 


and (B.2) 

V c i V c 2*z 

where ip ^ and \|> 2 are the atomic electronic wave functions of the projectile and 
the target, respectively. They are independent of the internuclear distance 
R. The coefficients C^, C^, C'p and C ' 2 are determined by the normalization 


condition and the LCAO method to be 
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Cl= (1 ♦ ESf^f 2 ) 1 ' 2 ; ^ 


C i ■ 


1 


(1 + 2Sfj+ f 2 ) 1/2 

frs 


o \ to ' » C 0 - 


(1 + 2Sf 2+ f^) 172 * 2 (1 + 2Sf 2+ f 2 ) 1/2 


where fj and f 2 are solutions of equations (2.18) and (2.20) and are known 
functions of the internuclear distance R. 

Now, first let us estimate how slowly the molecular eigen functions 
vary. For this consider 


V C lV C 2+2 


d C 1 dC„ 

= ^ dR *1 + ~~ dR~ * 7 1 - (£) 2 • 


where b is the impact parameter and is the velocity of the projectile. One 

notes that if is small, then ^ and similarly ^ should be small. Now, 
consider the following 


i = i A i|» e 
9 


- 1 “.Hi.."* J t£ S dt + Eg a* 9 e" 1 I*** dt 


+ i B t e 
u 


• 1 i' E ,« . t 8i e- 1 E* 


(B.3) 


\b 

ir u 


dt 


In the above, the second term contains a factor because of ^ ^ , and the 
third term contains a factor of E g which is about the order of unity. 
Therefore, if u 1, compared to the terms containing E g and E u , we can 
neglect the terms and . 
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Next, let us discuss the behavior of the quantities A(t) and B(t). For 
this we assume i » t u = 0, then using the Schroedinger equation it = Ht and 
equations (B.l) and (B.2), we obtain 

i A (C 1+ C 2 S) e-^ E g dt + A {Cj(E g - H n ) + C 2 (E g S - H 12 )} e _i J tE g dt 

+ iB (C{ + C’S) e _i ^ E u dt + B { C{(E U - H n ) + C‘ (E u S-H 12 )} e _i J tE u dt = 0 
Recall the basic equation of LCAO (eqs. (2.19) and (2.20)) 

C 1 (»!!” E ) + C 2^ H 12" eS ) = 0 


C 1 (H ?1 - e S) + C 2 (H 22 - O = 0 

and the two sets of solutions q, q and cq, Cq corresponding to eigen 
energies E g and E u , respectively. Thus, the equations 

Cj (Eg- H n ) ♦ C 2 (E g S - H 12 ) = 0 

C i <V V + C 2< E u S ' V ' 0 

lead to 

i A (Cj+ C 2 S) e" 1 ^ E g dt + i & (Cjq q S) e" 1 ^ E u dt = 0 (B.4) 
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Similarly, starting again from the Schroedinger equation = Hip, we will 
obtain 



The above two equations can be written as 

A (Cj+ C 2 S + & (CJ + C' 2 S) e' 1 V dt = 0 

and (8.5) 

A (C„+ C,S) + fc (Ci+ CIS) e“ 1 J'Vu- V dt = 0 

(B.6) 

Now, A and B have nonzero solutions if, and only if, the secular equation is 
equal to zero, as follows: 

(Cj+ C 2 S) (q+ C£S) e‘ 1 J l ( E u- E g } dt 

= 0 

(C 2 + CjS) (C’+ C'S) e” 1 JX- V dt 

which leads to 

V v (f 2~ f i )s2 

where S is the overlap integral and is always less than unity except when 
R = 0 and i|»j= <|> 2 , and where fj and f 2 , defined by fj = C 2 /C 1 and f 2 = C' 2 /C'p 
are obviously not always equal. Thus, the secular equation cannot always be 
zero. Therefore, the only solution of equations (B.5) and (B.6) is A = B = 0. 
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Thus, if we assume ij>g= u = 0, then A and B are independent of time. 

However, if u < 1, we can say that , A(t), and B(t) are very slowly 

9 u 

varying functions. In practice, if u < 5x 10 7 cm/sec , i. e., in atomic 
units u « 0.2 , the condition u < 1 is satisfied. It is noticed in figure 6 
for values u > 6 x 10 7 cm/sec, our results do not exhibit the increasing cross 
section as observed in experiment. 


TABLE 1. A list of functions for the reaction of 
Li + + H -► Li + (n=3) + H + 


FUNCTION 


EXPRESSION 


0.013 e‘ R (1.6 R 4 + 8R 3 + 24R 2 + 48R + 48) 

0.026 e _R (4R 3 + 8R 2 + 12R + 12) 

I - e- 2R ( 1 + I) 

R R 

- - (0.043 R 4 + 0.215 R 3 + 0.645R 2 + 1.29R + 1.61 + i) e' 

.R R 

- 0. 5S + — -K 


- .05 - 3.J + - 

R 

- 0.5 - J 2 + 1 


A i- (H 2 - H! ) 

Bi ( Jl + ( — ) 2 -1 1 A 

V 2A 

AE ( R ) (H 2 - Hj) + 2 Bj 


Hf 2 /(H 2 - Hi) 


(Hi*H 2 ) 


4 A 2 where, A = - ILiL£ + 

H 1 -H 2 -B 1 0.5 B 


zvr 


Hi = H 2 


V^+l 1 

(-) — 


(Hi- H 2 ) 
















Report Documentation Page 

bcace Ao^nisi'afor 


1. Report No. 2. Government Accession No. 

NASA TM-101559 


4. Title and Subtitle 

Charge Exchange Transition Probability for Collisions 
Between Unlike Ions and Atoms Within the Adiabatic 
Approximation 


3. Recipient's Catalog No. 


5. Report Date 

April 1989 

6. Performing Organization Code 


7. Author(s) 


Y. J. Xu, G . S. Khandelwal, J. W. Wilson 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665-5225 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


8. Performing Organization Report No. 

10. Work Unit No. 

199-22-76-02 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Y. J. Xu and G. S. Khandelwal: Old Dominion University, Norfolk, Virginia. 

J. W. Wilson: Langley Research Center, Hampton, Virginia. 


16. Abstract 

A simple formula for the transition probability for electron exchange between 
unlike ions and atoms is established within the adiabatic approximation by 
employing the Linear Combination of Atomic Orbitals (LCA0) method. The formula 
also involves an adiabatic parameter, introduced by Massey, and thus the 
difficulties arising from the internal en ^y^ y defect ^ yn d tljte adiab^yic appro^- 
mat|on are avoided. Specific reaction Li + H^Li + H and Be + H-^Be 
+ H are considered as examples. The calculated capture cross section results 
of the present work are compared with the experimental data and wjth the calcula- 
tion of other authors over the velocity range of 10 cm/sec to 10 cm/sec. 


17. Key Words (Suggested by Author(s)) 

Charge exchange 
Ion reactions 
Space Radiation 


18. Distribution Statement 

Unclassified-Unlimited 

Subject Category - 93 


19. Security Classrf. (of this report) 

20. Security Classif. (of this page) 

21. No. of pages 

22. Price 

Unclassified 

Unclassified 

42 

A03 


NASA FORM 1826 OCT 86 










